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We perform a global analysis of all available spin-dependent proton structure function data, 
covering a large range of Q 2 , 1 < Q 2 < 30 GeV 2 , and calculate the lowest moment of the <?i structure 
function as a function of Q 2 . From the Q 2 dependence of the lowest moment we extract matrix 
" elements of twist-4 operators, and determine the color electric and magnetic polarizabilities of the 

proton to be xe = 0.026 ± 0.015 (stat.) ± 8:824 (sys.) and \b = -0.013 =F 0.007 (stat.) =F 0.012 (sys.), 
^ ■ respectively. 
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Measurements of spin-dependent structure functions of the proton reveal fundamental information about the pro- 
ton's quark and gluon structure. In the quark-parton model, the gi structure function is interpreted in terms of 
distributions of quarks carrying light-cone momentum fraction x, with spins aligned versus anti-aligned with that of 
the nucleon. The lowest moment, or integral over x, of g\ also determines the total spin carried by quarks in the 
• nucleon. 

' Although most structure function studies have focused on the scaling regime at high four-momentum transfer 
Q\ , squared, Q 2 , the behavior of g\ and its moments in the transition region at intermediate Q 2 (~ 1 GeV 2 ) can reveal 
t-H ' rich information about the long-distance structure of the nucleon. One example of the complexity of this region is the 
transition from the Bjorken or Ellis-Jaffe sum rules at high Q 2 to the Gerasimov-Drell-Hearn sum rule at Q 2 = Q. 

Of particular importance in this region is the role of the nucleon resonances, and the interplay between resonant 
and scaling contributions. According to the operator product expansion (OPE) in QCD, the appearance of scaling 
violations at low Q 2 is related to the size of higher twist corrections to moments of structure functions Q • Higher twists 
' are expressed as matrix elements of operators involving nonperturbative interactions between quarks and gluons. The 
study of higher twist corrections thus gives us direct insight into the nature of long-range quark-gluon correlations. 

In this paper we determine the size of the higher twist contributions to the lowest moment of the g\ structure 
function of the proton. We analyze the entire set of available data from experiments at SLAC 0- EH- 0- E3l - CERN 0,0, 
1-^ 1 DESY and most recently Jefferson Lab where high-precision data in the resonance region and at low and 
intermediate Q 2 have been taken. Combining the moment data from the various experiments is nontrivial, however, 
since different analyses typically make use of different assumptions about extrapolations into unmeasured regions of 
kinematics. In the present analysis, we therefore extract the structure function moment using a single set of inputs 
and assumptions for all the data. 

The lowest (Cornwall-Norton) moment of the proton g\ structure function is defined as 

Ti(Q 2 )= f dx gi (x,Q 2 ) . (1) 
Jo 

The upper limit includes the proton elastic contribution at x = Q 2 /2Mv — 1, where v is the energy transfer, and 
M is the proton mass. The inclusion of the elastic component is essential if one wishes to use the OPE to study the 
evolution of the integral in the moderate Q 2 region . 

From the OPE, at large Q 2 the moment T% can be expanded in powers of 1/Q 2 , with the expansion coefficients 
related to nucleon matrix elements of operators of a definite twist (defined as the dimension minus the spin of the 
operator). At high Q 2 the moment is dominated by the leading twist contribution, [12, which is given in terms 
of matrix elements of the twist- 2 axial vector current, ip 7^75 ip. This can be decomposed into flavor singlet and 
nonsinglct contributions as 

^(Q 2 ) = c s (Q 2 )\ + c ns (g 2 )(g + I) (2) 

where C s and C ns are the singlet and nonsinglet Wilson coefficients, respectively 01 1 which are calculated as a series 
in a s . The triplet and octet axial charges, 0,3 = gA = 1.267 and as = 0.58, are extracted from weak decay matrix 
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elements. For the singlet axial charge, we work with the renormalization group invariant definition in the MS scheme, 
o,g nv = ao(Q 2 — oo), in which all of the Q 2 dependence is factorized into the Wilson coefficient C s . 

Considerable effort has been made over the past two decades in determining the singlet axial charge, which in the 
quark-parton model is identified with the total spin carried by quarks in the proton. In this work we focus instead 
on using the world's data to extract the coefficient of the l/Q 2 subleading, twist-4 term, which contains information 
on quark-gluon correlations in the nuclcon. 

In addition to the twist-4 matrix element, the l/Q 2 term also contains so-called "kinematic" higher twists, associated 
with target mass corrections (which are formally twist-2), and the g 2 structure function, which is obtained from 
measurements with transversely polarized targets. One technique for removing these from the l/Q 2 correction is to 
work in terms of the Nachtmann moment [Tflj . 

e J ( rf.fx 1 M 2 x£\ „ 2 ,4MV 



Mi(Q ) = / dx^t 9l (x,Q 2 )[---—^ -g 2 (x,Q 2 )-—^}, (3) 



where £ = 2x/(l+ y/l + 4M 2 x 2 /Q 2 ) is the Nachtmann scaling variable. The twist expansion of M\{Q 2 ) then yields 

M l( Q 2 ) = /, 2 (Q 2 ) + ^ + ^U--- (4) 

where /j,2 is given in Eq. J2J). The l/Q 2 correction in Eq. (0} exposes directly the "dynamical" twist-4 coefficient f%, 
since ^(Q 2 ) — 4f 2 (Q 2 ) /9M 2 , where f 2 is given in terms of a mixed quark-gluon operator, 

HQ 2 ) M 2 S» = i]Te 2 (N\g % G^ lv %\N). (5) 



Here G^ v — \^ va ^G a p is the dual of the gluon field strength tensor, is the proton spin vector, g is the strong 
coupling constant, and e q is the quark charge. 

Clearly the l/Q 2 term can be best determined in the intermediate Q 2 region, where Q 2 is neither so large as to 
completely suppress the higher twists, nor so small as to render the twist expansion unreliable. A meaningful analysis 
of data from different experiments further requires that the same set of inputs be used in the determination of gi, as 
well as g 2 . 

In practice one must reconstruct g\ from a combination of longitudinal (Ay) and transverse (A±) polarization 
asymmetries, together with the unpolarized F\ structure function, and the ratio R of the longitudinal to transverse 
cross sections. We begin by collecting all available data on An, as published in Refs. 0, 0> IE LJ IE HI IE US i and usc 
the same inputs in the analysis of A±, Fx and R for all the data sets. 

To provide a realistic description of Aj_, or equivalently A 2 (which is given in terms of Aii and A± 0), we 
consider both the resonance and nonresonant background contributions. For the background we use the (twist-2) 
Wandzura-Wilczek (WW) relation ^E]. Inclusion of target mass corrections in the WW formula 0] enables this 
prescription to be used down to low Q 2 , where target mass corrections are known to be important Q. In the 
resonance region, however, the WW approximation fails, and here alternative parameterizations are required. We 
calculate the resonance contribution from the electromagnetic helicity amplitudes Si/2(Q 2 ) and Ai/ 2 (Q 2 ) obtained in 
the constituent quark model [l5|. which includes 14 major resonances. The resonance contribution is then normalized 
to satisfy the Burkhardt-Cottingham sum rule fltl. The A 2 model is consistent with the available data [tH fl8j |. as 
well as with the model-independent Soffer bound |l9j . 

For the ratio R(x, Q 2 ) we use a new parametrization based on Rosenbluth-separated cross sections HE, El, which 
is adapted to the low-Q 2 and \ow-W 2 region, and smoothly interpolates to the earlier parameterization of the deep 
inelastic region from Ref. J^]. This parameterization uses all published data in the resonance region [2j|, as well as 
new data from Ref. [2p| . 

The spin-averaged proton cross section is well determined in both the resonance and deep inelastic regions. At high 
Q 2 the effect of R is small and the differential cross section da/dVldE' is proportional to F 2 (or Fi). At moderate Q 2 , 
however, the extraction of F\ and F 2 from the cross section requires knowledge of R. Using the parametrization of R 
from Ref. |2(} , we constructed a database of world data on F\ from which values of Fi corresponding to the measured 
Aii kinematic points were obtained by interpolation. Most of the data points for A\\ have kinematically close sets of 
d 2 a/dfldE' points, which allows interpolation uncertainties to be minimized. Full details of the extraction of g\ will 
be provided in a forthcoming publication [2l| . 

The resulting 171 structure function is illustrated in Fig. 1 as a function of x for several representative Q 2 values. 
The vertical arrows indicate the boundary between the resonance and deep inelastic regions sXW = 2 GeV. At the 
lower Q 2 values, Q 2 ~ 1 GeV 2 , a significant portion of the x range is in the resonance region (contributing ~ 40% of 
the magnitude of the lowest moment). 
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FIG. 1: Proton gi structure function at several different Q 2 values. The points represent reanalyzed experimental data obtained 
from the longitudinal asymmetry Ay from Refs. [il PI 1^11(1 FtIIM [ill 11^ using the procedure described in the text. The vertical 
arrows indicate the boundary between the resonance (to the right of the arrow) and deep inelastic regions (W = 2 GeV). 




FIG. 2: Q 2 dependence of the Nachtmann moment M\{Q 2 ). The error bars are statistical, with the systematic errors indicated 
by the hashed areas (see text). The leading twist (dashed), 1/Q 2 (dot-dashed), 1/Q 4 (dot-dot-dashed) and elastic (dotted) 
contributions are shown separately. The solid curve is the sum of leading and higher twist terms. 



The first moment of g\ is evaluated using the same method as in recent analyses of the unpolarized proton F2 
structure function moments pil I25I ] . This method is essentially independent of assumptions about the x dependence 
of the structure function when interpolating between data points, and is therefore well-suited for a study of Q 2 
evolution of the moments. For the low-x extrapolation, beyond the region where data exist, we use the Regge model- 
inspired parametrization from Ref. |26l. To estimate the uncertainty associated with the low- a; extrapolation, we also 
consider other parameterizations [27]], and take the maximum difference between the respective low- 2 contributions 
as the error. 

The resulting Nachtmann moment Mi(Q 2 ) is shown in Fig. 2, where the error bars on the data points are statistical 
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only. The systematic errors, some of which are correlated, are shown separately in the hashed areas above the data, 
and represent uncertainties from the low- a: extrapolation (lower hashed area) , and the experimental systematic errors 
together with those from A 2 , R and an estimated 5% uncertainty on the elastic contribution (upper hashed area). 
The g 2 contribution to M% is obtained from Am, A±, and Fi, as determined from the present analysis (see Ref. [21j 
for details). 

The fit to the total moment Mi(Q 2 ) uses three parameters, a™, f 2 (or /14) and fi6, with the non-singlet axial 
charges (gA and as) as inputs. For the leading twist contribution we use a next-to-leading order approximation for 
the Wilson coefficients and the two-loop expression for a s , which at Q 2 = 1 GeV 2 corresponds to a^ LO = 0.45 ± 0.05 
in the MS scheme. 

In fitting the parameters, we have considered both multiparameter (simultaneous) fits and sequential fits, in which 
the leading twist term a™ is first fitted to the high-Q 2 data, and then the higher twist terms are extracted. While 
both methods should in principle yield the same results when the experimental errors are small, in practice the multi- 
parameter fit may not be most suitable choice when emphasizing the high-precision low-Q 2 data. The multiparameter 
fit is most effective when the errors on the data are similar across the entire Q 2 range, and the number of points in 
the region which determines the leading twist contribution (Q 2 > 5 GeV 2 ) is similar to that which constrains the 
higher twists (Q 2 < 5 GeV 2 ). 

Assuming the data at high Q 2 are saturated by the twist-2 term alone, the fit to the Q 2 > 5 GeV 2 data determines 
the singlet axial charge to be 

afr y = 0.145 ± 0.018 (stat.) ± 0.103 (sys.) ± 0.041 (low x) ± 8;8?g (a s ) , (6) 

where the first and second errors are statistical and systematic, the third comes from the x — > extrapolation, and 
the last is due to the uncertainty in a s . We have considered the sensitivity of the results to the value of Q 2 used to 
constrain the leading twist term. We find that a™ v converges to the above value for Q 2 > 3-4 GeV 2 . Fitting the 
Q 2 > 10 GeV 2 data would lead to practically the same values of ag nv , but with a slightly larger error bar. 

Having determined the twist-2 term from the high-Q 2 data, we now extract the 1/Q 2 and 1/Q 4 coefficients from 
the 1 < Q 2 < 5 GeV 2 data, fixing aQ nv to the above value, but allowing it to vary within its statistical errors. For the 
twist-4 coefficient we find 

f 2 = 0.039 ± 0.022 (stat.) ± glSis MO ± °- 030 ( low x ) ± aoil K) . (<0 

normalized at a scale Q 2 = 1 GeV 2 (the Q 2 evolution of f 2 is implemented using the one-loop anomalous dimensions 
calculated in Refs. 28]). The systematic uncertainty on f 2 is determined by refitting the Mi data shifted up or 
down by the M\ systematic uncertainty shown in Fig. 2 (upper hashed area). The low- a: extrapolation uncertainty is 
determined by fitting the M\ values shifted by the maximum difference between the x — > contributions calculated 
with the parameterizations from Refs. [2(j and [27| (lower hashed area in Fig. 2). The relative contribution from the 
low- a; extrapolation to the asymptotic value of M\ is w 7% and 15% at Q 2 = 1 and 3 GeV 2 , respectively. However, 
because the x — > extrapolation affects more the overall magnitude of Mi rather than its Q 2 dependence, the effect 
on f 2 is relatively small. 

For the 1/Q 4 term, the best fit to the 1 < Q 2 < 5 GeV 2 data yields n 6 /M 4 = 0.011 ± 0.013 (stat.) ± g;ggg (sys.) ± 
0.011 (low a;)±0.000 (a s ), with the errors determined as for f 2 . Within the present level of accuracy the Q 2 evolution of 
the coefficient ^t 6 can be neglected. The 1/Q 2 and 1/Q 4 contributions to Mi are illustrated in Fig. 2. For comparison, 
the elastic contribution is also shown, as is the sum of the leading plus higher twist contributions. Since one is fitting 
Mi in a relatively low Q 2 region, one may ask whether still higher-order corrections could be significant beyond those 
considered in Eq. Adding a phenomenological /is/Q 6 term to the Q 2 > 1 GeV 2 fit, and fixing the other parameters 
to their quoted values, gives a coefficient Hs/M 6 = —0.004 ± 0.004 which is consistent with zero. To determine /xg 
more precisely one needs to go below Q 2 ~ 1 GeV 2 , however, it is not clear that the twist expansion is convergent at 
such low Q 2 . 

Simply fitting the entire 1 < Q 2 < 30 GeV 2 data set using a 3-parameter fit, the value of the singlet charge 
would be essentually unchanged (a™ = 0.145 ± 0.023 (stat.)), while the twist-4 coefficient would be slightly smaller, 
f 2 = 0.016 ± 0.039 (stat.), but compatible with the above result J7J within uncertainties (with similar systematic 
errors as in Eqs. © and (JJJ). 

In earlier phenomenological analyses [2^. [3 0ll larger values of f 2 were obtained. Usin g th e SLAC-E143 data 
Ref. [U found f 2 = 0.10 ± 0.05, while Ref. [3of used in addition HERMES and CLAS [l(J data to extract a value 
f 2 = 0.15 — 0.18. In both cases, however, the 1/Q 4 corrections in Eq. were not included, which we find important 
even for Q 2 ~ 1 GeV 2 . This is particularly relevant for the analysis in Ref. [3(j , which assumes that the higher 
twists are dominated by the 1/Q 2 term already at Q 2 ~ 0.5 GeV 2 . If one were to redo the present analysis with a 
1-parameter fit as in Refs. |29i l30| . the statistical error on f% from the Q 2 > 1 GeV 2 data would be ~ 0.005, which 
is 4-5 times smaller than that in the 2-parameter fit, and gives a 2-3 smaller combined statistical and systematic 



5 



uncertainty than in Ref. [2!j. The result of the present work J7J therefore represents a significant improvement over 
the earlier analyses. 

From the extracted f 2 values one can calculate the contribution of the collective color electric and magnetic fields 
to the spin of the proton. These are given by 0, 

XE = I (2rf 2 + h) , (8) 
Xb = \ (M 2 - h) , (9) 

where d 2 is given by the matrix element of the twist-3 operator i\;{G^ v ^ a + G fla 'y' / )ip, and can be extracted from the 
second moments of g\ and g 2 , 

d 2 (Q 2 ) = [ dx x 2 [2 gi (x,Q 2 ) + 3g 2 (x,Q 2 )], (10) 
Jo 

as determined from the present analysis. We find, however, that its leading twist component is negligible, and 
consistent with zero. 

Combining the extracted f 2 and d 2 values obtained from the global analysis, we find 

Xe = 0.026 ± 0.015 (stat.) ± (sys.) (11) 
Xb = -0.013 T 0.007 (stat.) T aoi" (sys.) (12) 

where the low-x extrapolation and a s uncertainties have been folded into the total systematic error. Since the color 
polarizabilitics are dominated by f 2 , the sign of the color electric polarizability is positive, while that of the color 
magnetic polarizability is negative. The upper limit on f 2 in Eq. thus yields non-zero values for \e and xb, while 
the lower limit gives values which are close to zero. 

These results can be compared to model calculations. QCD sum rules gen erally predict negative values for the 
electric polarizabilities, and slightly positive values for the magnetic ones [3l|, |33j, x^™ rulc ~ —(0.03 — 0.04) and 
Xb™ rule w 0.01 — 0.02, in contrast to the results in Eqs. (|llf> and (|12|) . Similar results are found in the calculations 
based on the instanton vacuum model 0, x'jf tanton w -0.03 and ^stanton « 0.015. The MIT bag model on the 
other hand gives 0] x^f 8 ~ 0.03 — 0.05 and x^ 8 ~ 0.00 — 0.02, which is consistent with our findings. 

More precise measurements of the structure functions at Q 2 « 1 — 10 GeV 2 , with smaller statistical and systematic 
errors, would reduce the uncertainty in the extracted higher twist coefficients, as would better knowledge of a s at the 
relatively low Q 2 values discussed here. The present findings suggest that higher twists in the lowest moment of the 
proton gi structure function are small and consistent with zero for Q 2 > 2 — 3 GeV 2 (see also Refs. H^,|3^), which 
demonstrates, perhaps surprisingly, the usefulness of the OPE formalism at these rather low Q 2 values. Higher twists 
are expected to play a greater role in higher moments, which emphasize more the large- x region and receive larger 
resonance contributions at the same Q 2 |2l| . Better determination of the g 2 structure function at moderate and high 
Q 2 is also vital for the determination of the d 2 matrix element, as well as of the g\ structure function itself. 
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DE-FG02-96ER41003 (College of William and Mary), and by NSF grant 0099540 (Hampton U.). 
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